The analytic structure of Bloch functions for linear molecular chains 
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This paper deals with Hamiltonians of the form H — — + v{r), with v{r) periodic along the 
z direction, v{x,y,z + b) = v{x,y,z). The wavefunctions of H are the weU known Bloch functions 
■!/;„, A(r), with the fundamental property ij)n,\{x,y,z + h) — \ipn,\{x,y, z) and dzilJn,\{x,y, z + 6) = 
^dz%l)n,\(x,y, z). We give the generic analytic structure (i.e. the Riemann surface) of 
and their corresponding energy, En{\), as functions of A. We show that En{X) and tp„^\{x,y, z) 
are different branches of two multi-valued analytic functions, E{X) and ^x{x, y, z), with an essential 
singularity at A = and additional branch points, which are generically of order 1 and 3, respectively. 
We show where these branch points come from, how they move when we change the potential 
and how to estimate their location. Based on these results, we give two applications: a compact 
expression of the Green's function and a discussion of the asymptotic behavior of the density matrix 
for insulating molecular chains. 

PACS numbers: 71.10.-w, 71.15.-m 



I. INTRODUCTION 

The analytic structure of the Bloch functions for ID 
crystals with inversion symmetry was investigated in 
Ref. Q. Among the major conclusions of this paper, is 
the fact that the entire band structure can be character- 
ized by a single, though multi-valued, analytic function 
E{X)[X = e^^^], with branch points that occur at com- 
plex k. A similar conclusion holds also for the Bloch 
functions. The positions of the branch points determine 
the exponential decay of the Wannier functions. For insu- 
lators, they also determine the exponential decay of the 
density matrix and other correlation functions. It was 
recently shown that the order of the branch points de- 
termines the additional inverse power law decay of these 
functions.'^ We can say that, although the branch points 
occur at complex their existence and location have very 
important implications on the properties and dynamics 
of the physical states. The complex k wavefunctions are 
also important when describing surface and defect states, 
metal-insulator junctions and electrical transport across 
finite crystals and linear molecular chainsi^i^*^ 

The methods developed in Ref. Q could not be ex- 
tended to higher dimensions, where the results are much 
more limited. The first major step here was made by 
Des Cloizeaux, who studied the analytic structure near 
real fc's, for crystals with a center of inversioni^ His 
conclusion was that the Bloch functions and energies of 
an isolated simple band are analytic (and periodic) in 
a complex strip around the real k's. The restriction to 
crystals with center of symmetry was later removedi^ The 
analytic structure has been reconsidered in a study by 
Avron and Simon^ who gave answers to some impor- 
tant, tough qualitative questions. For example, one of 
their conclusion was that all isolated singularities of the 
band energy are algebraic branch points. The topology 
of the Riemann surface of the Bloch functions for finite 
gap potentials in two dimensions has been investigated 
in a series of studies by Novikov et al.^^ The analytic 
structure has been also investigated by purely numerical 



methods. Since it is impossible to explore numerically the 
entire complex plane, numerical methods cannot provide 
the global structure. Even so, they can provide valuable 
information. For example, in the complex band calcu- 
lations for Si,^^ or linear molecular chains^i^ii^ one can 
clearly see how the bands are connected, even though 
these studies explored only the real axis of the complex 
energy plane. 

In this paper we consider linear molecular chains, de- 
scribed by a Hamiltonian of the form 



H 



v{v), 



(1) 



with v{y) periodic with respect to one of the cartesian 
coordinates of R'^, let us say z: 



v{x,y,z + b) = v{x,y,z). 



(2) 



The wavefunctions of H are Bloch waves, ipn^xir), with 
the fundamental property 

fpn,x{x,y,z + b) ^ Xijjn,\{x,y,z), (3) 
dz'il)n,x{x,y,z + b) = Xd^ipn,\{x,y,z). 

We derive the generic analytic structure (the Riemann 
surface) of i^n,\ (r) and of the corresponding energy En{X) 
as functions of complex A. We shall see that they are 
different branches of two analytic functions, and 
E(X), with an essential singularity at A = and addi- 
tional branch points which, generically, are of order 3, 
respectively 1. We show where this branch points come 
from, how they move when the potential is changed and, 
in some cases, how to estimate their location. 

Our strategy is as follows. According to Bloch theo- 
rem, we can restrict z to z G [0, b] and study the following 
class of Hamiltonians, which will be called Bloch Hamil- 
tonians: 

i^A = -V2 + i,(r), ze [0,&], (4) 
with A referring to the following boundary conditions: 



i'ix,y,b) = Xij}{x,y,Q), 
dzipix,y,b) = X^z^Pix,y,0), 



(5) 
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which define the domain of H\. For z € [0, 6], the eigen- 
vectors of Hx and their corresponding energies coincide 
with ^n,\{Y) and En{X). Let p{H\) denote the resolvent 
set of Hx, which is composed of those points in the com- 
plex energy plane for which {z^ H\)^^ is bounded. Now, 
it was long known that the Green's function, (z — H\)~^, 
evaluated at some arbitrary point z e p(Hx), is analytic 
of A, for any A in the complex plane»i2i In Section II, 
we derive the local analytic structure of the eigenvectors 
and eigenvalues as functions of A, starting from this ob- 
servation alone. To obtain the global analytic structure, 
we start with a simple v(r), with known global analytic 
structure, and then study how this structure changes 
when v{r) is modified. 

When this formalism is applied to ID periodic sys- 
tems, all the conclusions (and a few additional ones) of 
Ref. Q follow with no extra effort. This is done in Sec- 
tion III. Section IV presents the results for linear chains. 
We have two applications: a compact expression for the 
Green's function, which is developed in Section V, and 
the computation of the asymptotic behavior of the den- 
sity matrix for insulating linear molecular chains, which 
is done in Section VI. We conclude with remarks on how 
to calculate the analytic structure for real systems. We 
also have two Appendices with mathematical details. 



II. GENERAL FORMALISM 

We consider here, at a general and abstract level, an 
analytic family, {Hx}xgC, of closed, possibly non self- 
adjoint operators. The analyticity is considered in the 
sense of KatOjAfi. which means that, for any z E p{Hx), 
the Green's function can be expanded as 



{z-Hx'Y 



OO 



(A' - A)"i?„(z), 



(6) 



n=0 



with the power series converging in the topology induced 
by the operator norm, for any A' in a finite vicinity of A. 
We have already mentioned that the Bloch Hamiltonians 
form an analytic family. 

In this section, we discuss the analytic structure of 
the eigenvalues and the associated eigenvectors of Hx, 
as functions of A, based solely on the analyticity of the 
Green's function. For this, we need to find ways of ex- 
pressing the eigenvalues and eigenvectors using only the 
Green's function. The major challenge will be posed by 
the degeneracies. 

Suppose Hx has an isolated, non-degenerate eigenvalue 
Ex, for A near Aq. Then there exists a closed contour F 
separating Exg from the rest of the spectrum. For A in 
a sufficiently small vicinity of Ao, Ex remains the only 
eigenvalue inside F and we can express Ex as 



f , ^ 1 dz 

Ex^TrJ^ziz-Hxr-. 



(7) 




FIG. 1: The figure shows two isolated eigenvalues of H\, Ex 
and E'x, that become equal for A = Ac. The figure also shows 
the contours of integration, F, 7 and 7' used in the text. 



As shown in Appendix A, Eqs. © and O automatically 
imply that Ex is analytic at Aq. Since Aq was chosen ar- 
bitrarily, we can conclude that the non-degenerate eigen- 
values are analytic functions of A, as long as they stay 
isolated from the rest of the spectrum. 

Suppose now that Hx has two isolated, non-degenerate 
eigenvalues, Ex and E'^, which become equal at some Ac 
(see Fig. 1): 



Ex 



E' 



Er 



(8) 



We are interested in the analytic structure of these eigen- 
values and the associated eigenvectors, for A in a vicinity 
of Ac. 

Eq. Q is no longer useful, since there is no such A- 
independent contour F that isolates one eigenvalue from 
the rest of the spectrum, for all A in a vicinity of Ac. The 
key is to work with both eigenvalues, since we can still 
find a A-independent contour F, separating Ex and E'^ 
from the rest of the spectrum (see Fig. 1), as long as A 
stays in a sufRciently small vicinity of Ac. We define. 



F„,iX)^iExr' + iE'xr, m=l,2. 



(9) 



The main observation is that we can use the Green's func- 
tion to express F„i{X): 



Fm(A) =Tr / z"(z-i/A) 



dz 
27ri' 



(10) 



Then, as shown in Appendix A, it follows from Eqs. © 
and H10|l that F,„(A) are analytic functions near and at 
Ac. 

The functions introduced in Eq. (O provide the follow- 
ing representation: 



Ex 



E' 



F^[X) + ^2F2{\) - F^{\)A (11) 



Fi(A)- V2i^2(A)-Fi(A)^ 
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Thus, we managed to express the eigenvalues in terms of 
the Green's function alone. The analytic function, 

G(A) = 2^2(A)-Fi(A)2, (12) 

must have a zero at Ac, so its generic behavior near Ac is 

G(A) = (A-Ac)^.g(A), (13) 

with K an integer larger or equal to 1 and 5(A) analytic 
and non-zero at Ac- The cases when K > 2 are very spe- 
cial and will not be considered here. Only the following 
two possibilities are relevant to us: 

^ ' \ (A-Ac)^g(A) (type II), ^ ' 

where 5(A) is analytic and has no zeros in a vicinity of 
Ac. 

For a type I degeneracy, as A loops around Ac, E\ 
becomes E'^ and vice versa. The two eigenvalues are 
different branches of a double-valued analytic function, 
with a branch point of order 1 at Ac) 

Ex = E, + a{\ - Ac)i/2 ^ ^ ^ ^ ^ (15) 

For a type II degeneracy, both E\ and E'^ are analytic 
near Ac. 

We consieder now the spectral projectors. Pa and P^, 
associated with Ex and E'^, respectively. For A 7^ Ac, 
they have the following representation: 

Px = li^-Hx)-'^^, (16) 

where 7 is defined by |z — Ex\ — d, with d small enough 
(thus A dependent) so E'^ hes outside 7 (see Fig. 1). P'^ 
has a completely equivalent representation. We list the 
following properties: 

Pi = Px, P'x^ = P'x 
PxP'x = P'xPx = 0. (17) 

For a type I degeneracy, as A loops around Ac, Px be- 
comes P( and vice versa. Thus, Pa and Pi^ are differ- 
ent branches of a double-valued analytic function, with 
branch point of order 1 at Ac. In other words, they 
are given by the same function (identified from now on 
with Pa), which is evaluated on different Riemann sheets. 
These sheets have Ac as a common point. Thus, if Pa does 
not diverge at Ac, then 

hn. Pa = lirn P^. (18) 

A — *A(, A — *A(- 

But this will contradict, for example. Pa Pa = 0- We must 
conclude that Pa diverges at Ac. To find out the form of 
the singularity, we observe that, if /S.Ex denotes the dif- 
ference Ex — E'^, then AExPx has a well defined limit at 
Ac, which can be seen from the following representation: 

Ai?APA = l^{z - E'x){z Hx)-'^^. (19) 



Since A_Ea oc (A — Ac)^/^, we can conclude that the sin- 
gularity of Pa is of the form (A — Ac)^^/^. In addition, 
we mention that the total spectral projector, 

P(A)=PA + Pj., (20) 

is analytic near and at Ac, as it can be seen from the 
representation 

P{X)^ J^iz-Hxr'^^. (21) 

and that the Green function for A = Ac and z near Ec 
has the following structure: 

(z - HxJ-' - (z - E,)-^ lim AExPx 

+ {z-E,)-^P{K)+R{z), (22) 

with R{z) analytic near Ec- This expression is interesting 
because the Green's function of self-adjoint operators, 
viewed as functions of z, always have simple poles. Thus, 
Hx^ cannot be a self-adjoint operator. In other words, 
type I degeneracies cannot occur for those values of A for 
which Hx is self-adjoint. 

We now turn our attention to the eigenvectors. Since, 
for A 7^ Ac, Pa is rank one, it can be written as 

Px = \i^x){i^xl (23) 

with -ipx (V'a) the eigenvector to the left (right), normal- 
ized as 

(V^a,V'a) = 1. (24) 

As explained in Ref. passing from the projector to 
the eigenvectors is not a trivial matter, since these vec- 
tors are defined up to a phase factor. The question is, 
can we choose or define these phases so that no additional 
singularities are introduced? We can give a positive an- 
swer when there is an anti- unitary transformation, Q, 
such that: 

H{ = QHxQ-\ (25) 

As we shall later see, such Q exists, for example, when 
v{x,y,z) — v{x,y, —z). Now fix an arbitrary ip and ob- 
serve that 

Pa ^ E^}fimE^. (26) 

(p|q^,PaV> 

One immediate problem with the above expression is that 
the denominator can be zero. This can happen only for 
isolated values of A for, otherwise, the denominator will 
be identically zero. Let Ao be such value, assumed dif- 
ferent from the branch point (we can always choose a ip 
satisfying this condition). We know that the only singu- 
larity of Pa is at Ac, so what we have is two functions that 
are equal on a domain surrounding Ag but excluding Aq, 
and one of them. Pa, is analytic at Aq. Then it is a fact 
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that both functions are analytic at Aq. This means the 
numerator in the right hand side of Eq. I|26|) must also 
cancel at Aq, with the same power as the denominator 
and the problem disappears. Then, it is natural to think 
that we can define the left and right eigenvectors of Hx 
as: 



and 



(27) 



(28) 



Using the properties of Q, we can equivalently write: 



I^a) = Q|^a). (29) 



The only problem is the following: the denominator in 
Eq. H26|l behaves like (A — Aq)" near Aq, with a some 
unknown power. As we already seen, the numerator of 
Eq. H26|l must have exactly the same behavior. Now let 
us look at Eq. (|29l) : the denominator of tpx behaves like 
(A — Ao)"/^, but the numerator can behave like (A — Ao)'^ 
with r arbitrary, as long as the numerator of ipx behaves 
like (A — Ao)""'^. However, we show in the following the 
r is exactly a/2. If Q cannot be defined, then this is 
no longer true. You may think that, in such cases, we 
should modify the denominator of ^px to {Q Pxip , Px'^'V ■ 
The problem is that Ao is not unique so this may fix the 
problem at Aq but not at other similar points. 
Now let 



(30) 



be the expansion of Px near Aq. If the numerator of 
Eq. (|26|l cancel at Aq, an expansion in powers of A will 
show that this cancelation is equivalent to: 

0, for alH < /s:, (31) 

with K an integer larger or equal to 1. This also means 

oo 

^aI^) = 5](A-Ao)'P.|V) = (A-Ao)^|0a). (32) 

i=K 

Let us assume, for the beginning, that K = 1, in which 
case Pi IV') 7^ 0. The numerator of Eq. H2t)|) becomes 

\PxtP){Pxm = {X-Xof\Piij){Piij\Q + o{X~Xof (33) 

and the denominator 

{QPxi^,Px'4>) = {X-Xof{QPii',Pii')+oiX-Xof. (34) 

Although 7^ 0, we are not automatically guaran- 

teed that {QPi4>,Piip) 0. However, we already argued 



that numerator and denominator must cancel with the 
same power, so this must be so. We can repeat the same 
arguments for arbitrary K and the conclusion will be the 
same: 



{QPx^,Pxi^) = 



(A- 



Ao)'^5(A) 
-Ao)^|0a), 



(35) 



with g{X) and \(j)x) analytic and non-zero near Aq. Thus, 
we can take the square root in Eq. (|29|) without intro- 
ducing a branch point. We can also see that the denom- 
inator and numerator in Eq. (|29|) cancel at Aq with the 
same power, so there is no pole at Aq. The conclusion is 
that Ipx is analytic at Aq. 

There will be, inherently, a branch point at Ac, where 
V'A and ■ipx behave as 



V'A(r) 
rxir) 



(A-Ac)V4 

5(r) 

(A - Ac)i/4 



■d(r) 
■J(r) 



(36) 
(37) 



i.e. ipx and ipx have a branch point of order 3 at Ac. If the 
operator Q with the above mentioned properties exists, 
this is their only singularity near Ac. 

For a type II degeneracy, Px and P'^ are analytic near 
Ac. The eigenvectors can be introduce in the same way 
as above and, if Q exists, they are analytic functions of 
A. 

Analytic deformations. We analyze now what happens 
when an analytic potential w is added: 



Hx,j = Hx+ -fw. 



(38) 



By analytic potential we mean that {Hx.'y}x,'yec is an 
analytic family in the sense of Kato, in both A and 7 (see 
Appendix B). 

For any fixed 7, the isolated, non-degenerate eigenval- 
ues of ffA,7 remain analytic of A. The interesting question 
is what happens with the degeneracies. Suppose that, 
at some fixed 70, there are two isolated, non-degenerate 
eigenvalues, -Ea,7o and E'^ which become equal at Ac. 
For A in a small vicinity of Ac and 7 in a small vicinity 
of 7o, we can define the functions Fi^2(A,7) and G(A,7) 
as before, which are now analytic functions in both ar- 
guments, (A, 7), near (Ac,7o). 

If at 7o , Ac is a type I degeneracy, the only effect of a 
variation in 7 is a shift of Ac. Indeed, since 



G(Ac,7o) = 0, aAG(Ac,7o) y^O, 



(39) 



the analytic implicit function theorem assures us that 
there is a unique Ac (7) such that 



G(Ac(7),7) = 0. 



(40) 



Moreover, the zero is simple, i.e. Ac(7) remains a type I 
degeneracy. For 7 near 70, 

Ac(7)^Ac- (7-7o)-K..., (41) 

CAG(Ac,7oj 



where the dots indicate higher order terms in 7 — 70- 
From Eqs. (|15|l and l|3t)|) we readily find: 



Ac(7) = A, 



I _ 2. — 72. j c{x)w{x)c{x)dx + 



(42) 



If iJA,7 and -ffA*,7 have complex conjugate eigenvalues, 
then: 



G(A,7) = G(A*,7)* 



(43) 



In this case, if Ac is located on the real axis, so it is Ac (7). 
If not, then 



G(Ac(7),7)-G(Ac(7)*,7) = 0, 



(44) 



which contradicts the uniqueness of Ac (7). 

Generically, a type II degeneracy splits into a pair of 
type I degeneracies when 7 is varied. Indeed, since 

G(Ac, 70) = dxG{\c, 70) - a^G(Ac, 70) = 0, (45) 

the generic structure of G(A,7) near (Ac, 70) is: 

G(A,7) =a(A-Ac)2+6(A-Ac)(7-7o) + c(7-7o)' + .... 

(46) 

Thus, for 7 7^ 7o, G(A,7) will have, generically, two sim- 
ple zeros (type I degeneracies) at: 



A?(7) = Ac + 



7-70 
2a 



+ .. 



(47) 



The coefficients a, b and c can be derived from a pertur- 
bation expansion of Eqs. lO and H12() . leading to 



Tr {(Pa. - P'xM ± J -Tr{Px^wP;,w} , (48) 



plus higher orders in 7 — 70. 

If a type II degeneracy does not split, it remains a type 
II degeneracy for all values of 7. This is a consequence 
of the fact that, if two analytic functions are equal on an 
interval, they are equal on their entire domain. 

We now summarize the findings of this Section. The 
isolated, non-degenerate eigenvalues of H\ are analytic of 
A. Double degeneracies can be of different kinds. Type I 
degeneracies are equivalent to branch points. Near such 
degeneracies, E\ behaves as a square root. Type I de- 
generacies are robust to analytic perturbations: as long 
as they stay isolated, variations of the periodic poten- 
tial cannot destroy or create but only shift them. At 
type II degeneracies, the eigenvalues are analytic. Type 
II degeneracies are unstable to analytic perturbations: 
generically, they split into two type I degeneracies. The 
other types of degeneracies were considered rare and not 
discussed here. The analytic structure of the spectral 
projectors can be automatically deduced from the ana- 
lytic structure of E(X). If there exists an anti-unitary 
transformation, Q, such that h\ = QH\Q~^, then the 





(b) 



FIG. 2: a) The Riemann surface of E^{\) is a spiral, which 
have been cut along the dotted lines in individual sheets, in- 
dexed by n = 0, ±1, . . .. The solid dots indicate the type II 
degeneracies and the arrows indicate how they pair, b) The 
Riemann surface of E{\) at 7 > 0. The empty dots represent 
the branch points and the arrows indicate how the Riemann 
sheets are connected. In both panels, the thick line shows the 
trajectory on the Riemann surface, when A moves on the unit 
circle. 



phase of the eigenvectors can be chosen in a canonical 
way and their analytic structure follows automatically 
from the analytic structure of E{\). If such Q does not 
exists, the eigenvectors will still have singularities of the 
type (A — Ac)~^/^ at the branch points of E\ but the 
present analysis does not rule the existence of additional 
singularities. 

A similar analysis can be developed for higher degen- 
eracies. For a triple degeneracy, for example, we will 
have to consider the functions Fm{X), with m = 1,2,3. 
However, we regard higher degeneracies as non-generic, 
i.e. more rare than simple degeneracies and will not be 
considered in this paper. 



III. STRICTLY ID SYSTEMS 

We apply here the abstract formalism to an already 
well studied problem: the analytic structure of Bloch 
functions in ID, i.e the wave functions of the following 
Hamiltonian: 



H = —d^ + v{x), v{x + b) = v{x), X £ 



(49) 



According to Bloch theorem, finding the wave functions 
and their corresponding energies is equivalent to studying 
the following analytic family of Hamiltonians:^^ 



Hx = -{dl)^+v{x), xe[0,b], 



(50) 



defined in the Hilbert space of square integrable functions 
over the interval [0, b] . The index A refers to the boundary 



FIG. 3: a) The plot of ^^"(A) (n = 0, ±1, . . .) as functions of 
real fez (A = e'*^'). b) The generic band structure when the 
periodic potential is turned on. 

conditions 

m = xm, i^'ip) = xi^'iQ), (51) 

which define the domain of H\. The energy spectrum 
of H consists of all eigenvalues of Hx^ when A sweeps 
continuously the unit circle. li^n,\{x) is the normalized 
eigenvector of H\ corresponding to the eigenvalue i?„(A), 
then ipn,\{x) coincides on [0, b] with the Bloch wave of the 
same energy. If we extend these functions to the entire 
real axis by using 

^„,a(x + mb) = A'"V„,a(x), X e [0, b], (52) 

they will automatically satisfy the standard normaliza- 
tion, 

/oo 
'4>n,i/\{x)il^n,\'{x)dx = 2'KiXS{X - A'). (53) 
-oo 

Here are a few elementary properties of Hx- H\ is an 
analytic family in the sense of Kato, for all A G C. Hx 
is self-adjoint if and only if A is on the unit circle. In 
general, -ffi/A* is the adjoint of Hx- If C is the com- 
plex conjugation, CHxC — Hx*. Thus, Hx and Hx" 
have complex conjugated eigenvalues and Hx and Hi/x 
have identical eigenvalues. This tells us that the analytic 
structure is invariant to A — > 1/A and A — > A*. Because 
of these symmetries, it is sufficient to consider only the 
domain |A| < 1. For A not necessarily on the unit circle, 
the spectral projector on E„ x is given by 

PnA^' y) = ^n,x{x)'4'n,l/\{y)- (54) 

We consider now the following class of Hamiltonians 

Hxrt^-{dl)^+-iv{x), v{x + b):^v{x), (55) 

and we adiabatically switch 7 from zero to one. As it 
is shown in Appendix B, potentials v{x) with square in- 
tegrable singularitiesi^ are analytic. Thus, the theory 



FIG. 4: An equivalent representation of the Riemann surface 
of Fig. 2b. Each sheet corresponds now to one energy band. 
The thick line shows the trajectory on the Riemann surface 
when A encircles the origin at a small radius. 

developed in the previous Section can be applied to a 
large class of potentials. 

The eigenvalues and the associated eigenvectors of 
Hxfi are given by: 

El{X) = 6-2(2n7^i + lnA)^ 

V,0,^(x) = 5-l/2g(2«.«+l„A)./6_ (5g^ 

They are the different branches of the following multi- 
valued analytic functions: 

E"{X) = (inA)2 , ^jI{x) = 6-1/2 A"/ft. (57) 

The Riemann surface of E'^{X) is shown in Fig. 2a. There 
are only type II degeneracies and they occur at A — ±1: 

i?°(A = l) = i?^(A = l), 

E^,{X = -1) = i?Vi(^ = -l)- (58) 

A plot of these eigenvalues for A on the unit circle is 
given in Fig. 3a. Now we turn on the periodic potential. 
Typically, the energy spectrum of H consists of an infinite 
set of bands separated by gaps, denoted here by A„ (see 
Fig. 3b). When v{x) is modified, some of the gaps may 
close and other may open. If we assume that all the 
gaps open when the periodic potential is turned on, as 
in Ref. nj, then all type II degeneracies split into pairs of 
type I degeneracies, Ac (7) and 1/Ac(7), located on the 
real axis. Ac (7) and 1/Ac(7) branch points for E{X), 
which connect different sheets of the original Riemann 
surface (see Fig. 2b). Since they are constrained on the 
real axis, the trajectory of different branch points cannot 
intersect (they stay on the same Riemann sheet), i.e. the 
branch points remain isolated as 7 is increased. Thus, as 
the previous Section showed, they move analytically as 
we increase the coupling constant and we can conclude 
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\m(E-E.)\^i(E-E.) 




FIG. 5: /i(£' — Ej) and ^{E — Ei) as functions of E, for a 
typical Kramers function /i. The solutions to Eq. 1691 are 
given by the intersection points A, B, C, .... 



that the analytic structure cannot change, quahtatively, 
as we increase 7. 

For 7 = 0, we move from one sheet to another as A 
moves continuously on the unit circle, as Fig. 2a shows. 
The situation is different for 7 7^ (see Fig. 2b): when A 
completes one loop on the unit circle, we end up on the 
same point of the Riemann surface as where we started. 
We can then re-cut the Riemann surface, so that we stay 
on the same sheet when A moves on the unit circle (see 
Fig. 4). 

Thus, we rediscovered one of the main conclusions of 
Ref. Q. The eigenvalues En{\) are different branches of 
a multi-valued analytic function E{X) with a Riemann 
surface shown in Fig. 4: there are branch points or order 
1 at Ai, A2, . . ., and an essential singularity at 0. For 7 
small, Eq. (gSJl leads to 



A„ = (-1)'^ 



1 - 



6A„ 



(59) 



where A„ is the n-th energy gap and e„ is the energy in 
the middle of the gap. 

If v(x) = v{—x), we can construct Q as Q — CS, 
where C is the complex conjugation and S is the inver- 
sion relative to x = 0. Thus, for systems with inver- 
sion symmetry, we can also conclude at once that the 
only branch points of the Bloch functions are Ai , A2 , . . . , 
which are of order 3 (see Eqs. (j^ ). The present analy- 
sis actually adds something new to the results of Ref. 1., 
where the author studied two particular phase choices 
of the Bloch functions, namely, those imposed by ip\{x) 
(|A| ~ 1) being real at the points of inversion symme- 
try, X — and x — 6/2. The author warns that other 
choices can introduce additional singularities and thus 
reduce the exponential localization of the corresponding 
Wannier functions. A frequently used method of gener- 
ating Wannier functions localized near an arbitrary xq is 
to impose lmip\{xo) = 0. Since such a phase choice cor- 
responds to choosing ip{x) = d{x — Xo) in Eq. H29|l . we can 




FIG. 6: a) Ei_. (A) and E^. (A) as functions of real (A = 
gifczt^ for case A. The dashed lines shows the bands after 
the non-separable potential was turned on. b) The Riemann 
sheets of . (A) and i?^ . (A) and the type II degeneracies 
(solid circles), with arrows indicating how they pair (case A), 
c) The Riemann surface at 7 > 0. The empty circles represent 
the branch points and the arrow indicate how they connect 
different points of the Riemann surface. 



automatically conclude that it does not introduce addi- 
tional singularities and that the exponential localization 
of the corresponding Wannier functions is maximal. 



IV. LINEAR MOLECULAR CHAINS 

We specialize our discussion to 3D and consider Hamil- 
tonians of the form: 

H = -V^ + V{r), Vix,y,z + b) = Vix,y,z). (60) 

Using the Bloch theorem, we can find the spectrum and 
the wave functions of H by studying the following family 
of analytic Hamiltonians: 



^dl-dl-id,)l + V{r), [0,5], (61) 



with the boundary conditions 



i^{x,y,b) = AV'(a;,y,0) 
d^il}(x,y,b) = Xd:,ip{x,y,0). 



(62) 
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FIG. 7: The Riemann sheets of -B±(A) (case A). The thick hne 
shows a trajectory on the Riemann surface when A completes 
a loop around the origin. 

The general facts about Hx listed in the previous Section 
are still valid. Again, the analytic structure is invariant 
to A ^ A* and A so we can and shall restrict the 

domain of A to the unit disk, |A| < 1. 

We apply the analytic deformation strategy, as we did 
for the ID case. We start from a Hamiltonian with known 
global analytic structure. For this we consider a separa- 
ble potential, 

Hn = -V^ + v^{x, y) + v{z), v{z + b) = v{z), (63) 

and then adiabatically introduce the non-separable part 
of the Hamiltonian, 

= Ho+"fw{v), w{x,y,z + b) ^w{x,y,z). (64) 

We assume, for simplicity, that 

H^ = -dl-dl+VA^ix,y) (65) 

has only discrete, non-degenerate spectrum. For this, we 
will have to constrain x and 1/ in a finite region, which 
can be arbitrarily large. To be specific, we assume x,y d 
[0, 6'] and impose periodic boundary conditions. This 
is actually the most widely used approach in numerical 
calculations involving linear chains. We also assume that 
w(r) has square integrable singularities, which warranties 
that it is an analytic potential (see Appendix B). 

Let 4ij{x,y), Ej and ipn^\{z), £'„(A) denote the eigen- 
vectors and the corresponding eigenvalues of H± and of 

H\\x^-{d,)l+v{z), (66) 

respectively. Then the eigenvectors and the correspond- 
ing eigenvalues of Hq are given by: 



(a) 




FIG. 8: a) E^^{X) and i?^. (A) as a function of real 
(A = e^^'^) for case A. The dashed lines shows the bands 
after the non-separable potential was turned on. b) The Rie- 
mann sheets of i?^ . (A) and i?^. (A) and the type II degenera- 
cies, with arrows indicating how they pair (case C). c) The 
Riemann surface at 7 > 0. The empty circles represent the 
branch points and the arrow indicate how they connect dif- 
ferent points of the Riemann surface. 

The global analytic structure of i?^(A) is known: for j 
fixed, they are different branches of a multi-valued ana- 
lytic function E^ {X), with a Riemann surface as in Fig. 4. 
We now look for type II degeneracies: 

E = E,+ En^ {X) = E,+ En^ (A) , (68) 

which can occur only for A on the unit circle or on the 
real axis but away from the branch cuts. Indeed, if /i(i?) 
denotes the Kramers function for the strictly one dimen- 
sional Hamiltonian H^^xr^ then Eq. H68|l is equivalent to 

^iiE-Ej)=fi{E-Ei). (69) 

In Ref. ilj, it was shown that the equation dfi/dE = has 
solutions only for E on the real axis. Using exactly the 
same arguments, one can show that all the solutions of 
Eq. H69|) are on the real axis (see Fig. 5). Given that 

X^ -2^i{E-Ej)X + l = Q, (70) 

with fx{E — Ej) real, it follows that A must lie either 
on the unit circle or on the real axis (away from the 
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FIG. 9: The Riemann sheets of E±{\) (case C). The thick 
hnes shows a trajectory on the Riemann surface when A com- 
pletes a loop around the origin. 

branch cuts). For example, the solutions A and C in 
Fig. 5 have A on the unit circle, while the solution B 
has A on the real axis, inside the unit circle. We will 
refer to this three situations as cases A, B and C. Since 
the analytic structure is symmetric to A ^ A*, the type 
II degeneracies on the unit circle always come in pair, 
symmetric to the real axis. 

We consider first the case A, which corresponds to the 
case when two bands intersect as in Fig. 6a. Fig. 6b 
shows the Riemann sheets corresponding to these two 
bands. There are two type II degeneracies, marked with 
solid circles, on the unit circle and symmetric to the real 
axis. We assume for the beginning that these are the only 
type II degeneracies that split when the non-separable 
potential is turned on. When the type II degeneracies 
split, avoided crossings occur and the bands split (see 
the dashed lines in Fig. 6a). We denote the upper/lower 
band by E±{X). When the non-separable part of the po- 
tential is turned on, the already existing branch points 
shift along the real axis. For small 7, the shifts can be 
calculated from Eq. H42() . In addition, two type I degen- 
eracies appear (and another two outside the unit disk), 
introducing branch points that connect the original Rie- 
mann sheets. The connected sheets are shown in Fig. 6c, 
where we can also see that, when A makes a complete 
loop on the unit circle, we end up at the same point as 
where we started. This means we can re-cut the two, 
now connected, sheets so that we stay on the same sheet 
when A moves on the unit circle. These new sheets are 
shown in Fig. 7 and correspond now to the upper/lower 
bands E±{X). 

We consider now the case C, which corresponds to a 
situation when two bands intersect as in Fig. 8a. When 
the type II degeneracies split, the bands split in E±{X) 
and a gap appears. This is the only qualitative differ- 
ence between cases A and C. The Riemann surfaces, be- 
fore and after the non-separable potential was turned on. 




FIG. 10: a) The band Ei^{X) intersects with K^(A) and 
E^^{X). The dashed lines shows the bands at 7 > 0. b) The 
Riemann sheets for these bands and the type II degeneracies, 
with arrows indicating how they pair, c) The Riemann sheets 
at 7 > 0. 

are shown in Figs. 8b and 8c. Again, we can re-cut the 
Riemann surface so one sheet corresponds to one band. 
These new Riemann sheets are shown in Fig. 9. 

The case B goes completely analogous. The qualitative 
differences are that the branch points split from the real 
axis and we don't have to re-cut the Riemann surface. 

We analyze now a more involved possibility, namely 
when we have more type II degeneracies on the same 
Riemann sheet: 

E^X)=e:^{X), E^X')^Eliy). (71) 

Such a situation appears when, for example, we have 
bands crossing as in Fig. 10a. The Riemann sheets for 
these bands and the type II degeneracies are shown in 
Fig. 10b. When the perturbation is turned on, the type 
II degeneracies split in pairs of type I degeneracies, in- 
troducing branch points connecting the original sheets as 
shown in Fig. 10c. Again, when A completes one loop on 
the unit circle, we end up at the same point of the Rie- 
mann surface as where we started. We can then re-cut 
the Riemann surface such that we stay on the same sheet 
when A moves on the unit circle (see Fig. 11). The only 
new element is a Riemann sheet (corresponding to the 
middle band) with 6 branch points. 
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FIG. 11: The Riemann sheets corresponding to the three 
bands (dashed hnes) of Fig. 10a, for 7 > 0. 



The last situation we consider is the emergence of a 
complex band. Suppose that H has a symmetry with an 
irreducible representation of dimension 2. Suppose that 
this symmetry is also present for the Bloch Hamiltonian 
at A = 1. In this case, the separable Hamiltonian will 
have bands that touch like in Fig. 12a. Such situations 
are no longer accidental. In this case, the function G(A) 
introduced in Eq. H12|l behaves as 



G(A) = (A-1)VA), 



(72) 



with g{X) non-zero at A = 1. Following our previous no- 
tation, this will be a type IV degeneracy (see Fig. 12b). 
When the non-separable potential is turned on, the de- 
generacy at A = 1 cannot be lifted because of the sym- 
metry. This means G(A) must continue to have a zero at 
A = 1. Generically, the order of this zero is reduced to 2 
and two other zero's split, symmetric relative to the unit 
circle. In other words, the type IV degeneracy splits into 
a pair of type I degeneracies plus a type II degeneracy. 
E{X) remains analytic at A = 0, but now the two bands 
are entangled, in the sense that we need to loop twice on 
the unit circle to return back to the same point of the 
Riemann surface (see Fig. 12c) and we can no longer cut 
the Riemann surface so that we stay on the same sheet 
when A moves on the unit circle. A complex band can 
involve an arbitrary number of bands. For example, the 
new bands shown with dashed lines in Fig. 12a can en- 
tangle with other bands at A = —1, through the same 
mechanism, and so on. Rather than cutting the Rie- 
mann surface in individual sheets, we think it is much 
more convenient to think of a complex band as living on 
a surface made of all the individual sheets that are en- 
tangled through the mechanism described in Fig. 12. We 
note that the complex band can split in simple bands as 
soon as the symmetry is broken. 

We can continue with further examples but we can al- 
ready draw our main conclusions. The eigenvalues of 
i?A,7 are different branches of a multi-valued analytic 
function E{X). The Riemann surface of £'(A) can be cut 



in sub-surfaces, such that each sub-surface describes one 
band. For a simple band, this subsurface consists of the 
entire unit disk, with cuts obtained by connecting a fi- 
nite number of branch points to the essential singularity 
at A = 0. For a complex band, the sub-surface consists 
of a finite number of unit disks that are connected as in 
Fig. 12c. On this sub-surface we can have an arbitrary 
number of branch points, that connect this sub-surface 
to the rest of the Riemann surface. For both simple and 
complex bands, the branch points are symmetric relative 
to the real axis and, generically, they are of order 1 ("ac- 
cidental" higher degeneracies can lead to branch points 
of higher order). 

As we analytically deform the Hamiltonian, the un- 
split type II degeneracies stay on the unit circle or real 
axis and the positions of the branch points shift smoothly. 
In contradistinction to the ID case, the branch points can 
move from one sheet to another and their trajectories 
can intersect. When two of them intersect, they either 
become branch points of order 2 or recombine into a type 
II degeneracy. Higher order branch points are not stable, 
in the sense that small perturbations split them into two 
or more branch points of order 1. 

The Riemann surface of the spectral projector Px is 
the same as for E{X). When the inversion symmetry is 
present, the analytic structure of the Bloch functions can 
also be completely determined from the analytic struc- 
ture of E{\): the Riemann sheets of f^re the same 
as for E{X), but the branch points are generically of order 
3 (see Eq. (ED)). 



V. THE GREEN'S FUNCTION 

With the analytic structure at hand, it is a simple exer- 
cise to find a compact expression for the Green's function 
Ge = [E — H)^^, which is a generalization of the well 
known Sturm-Liouville formula in ID. Indeed, using the 
eigenfunction expansion. 



GB(r,r') = 



V'„4/A(r)V'n,A(r') dX 



^ '|A| = 1 ^ - En(X) 



27riA' 



(73) 



where the sum goes over all unit disks of the Riemann 
surface. Changing the variable from A to 1/ A if necessary, 
the above expression can also be written as: 



V-'«4/A(r<)'0n,A(r>) dX 



|A| = 1 



EniX) 



2mX' 



(74) 



where r> = r if 2: > z' , r> = r' if z' > z, and sim- 
ilarly for r< . This step is necessary because we will 
deform the contour of integration inside the unit circle. 
We could deform the contour outside the unit circle, but 
since we chose to exclude this part of the domain we 
don't have this liberty anymore, and we need to re-arange 
the arguments before deforming the contour. The inte- 
grand (including the summation over n) is analytic at 
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and un-occupied states. We start from 




FIG. 12: a) Two bands, E'„.{X) and Ei.{\), touch tangen- 
tially at fcz = (A = 1). The dashed hnes show the bands 
after the non-separable potential was turned on. b) The Rie- 
mann sheets for these bands and the type IV degeneracy (solid 
circle), c) The Riemann sheets after the non-separable poten- 
tial was turned on. 



the branch points. Also, for A ^ 0, En{X) — > oo and 

'0n,i/A(r<)V'ii,A(r>) oc A'^^^ I, so there is no singularity 
at A = 0. Then, apart from poles, which occur whenever 
En{X) — E, the integrand is analytic. Using the residue 
theorem, we conclude 



V'i/A,(r<)V'A, (r>) 
XjdxEiX,) ■' 



(75) 



where the sum goes over all Xj on the Riemann surface 
such that E{Xj) = E. This expression is valid for sys- 
tems with and without inversion symmetry, since it is the 
projector, not the individual Bloch functions, that enters 
into the above equations. Eq. (|75|l is closely related to the 
surface adapted expression of the bulk Green's function 
derived in Ref. 1^ 



VI. THE DENSITY MATRIX 

As a simple application, we derive the asymptotic be- 
havior of the density matrix n(r,r') for large \z — z'\, 
when there is an insulating gap between the occupied 



n(r, r') 



/ GE{r,r')dE, 
Im Jr 



(76) 



where C is a contour in the complex energy plane sur- 
rounding the energies of the occupied states. Using 
Eq. H75|l . we readily obtain 



'^(r,r') = ^ y" V'i/A(r<)V'A(r>)A ^dA, 



(77) 



where 7 is the pre-image of the contour C on the Riemann 
surface of E{X). We now restrict r and r' to the first unit 
cell and calculate the asymptotic form of n(r, r' + mbez) 
for large m. Using the fundamental property of the Bloch 
functions, we have 

n(r,r' + m6e,) = ^ ^i/;,(r<)^;,(r>)A"-idA. (78) 

We deform the contour 7 on the Riemann surface such 
that the distance from its points to the unit circle is max- 
imum. In this way, we enforced the fastest decay, with 
respect to m, of the integrand. This optimal contour, 
will surround (infinitely tide) the branch cuts enclosed 
by the original contour. The asymptotic behavior comes 
from the vicinity of the branch points Ac and A* (they 
always come in pair) that are the closest to the unit cir- 
cle. Using the behavior of the Bloch functions near the 
branch points, we find 



Re c(r)c(r')A^ 



(A/Ac)'"-M(A/Ac) 
^(l-A/Ac)i/2 ' 



(79) 



where the integral is taken along the branch cut of Ac. 
This integral is equal to ^B{m, 1/2), with B the Beta 
function. Wc conclude: 

n(r, r' + mbe,) -B{m, l/2)Re[c(r)c(r')A^]. (80) 

TT 

Again, this expression holds for systems with and with- 
out inversion symmetry, since it is the projector, not the 
individual Bloch functions, that enters in the above equa- 
tions. 



VII. CONCLUSIONS 

First, we want to point out that the formalism pre- 
sented here can be also applied to cubic crystals, to de- 
rive the analytic structure of the Bloch functions with 
respect to fcz, while keeping kx and ky fixed. Prelimi- 
nary results and several applications have been already 
reported in Ref. |^ 

We come now to the question of how to locate the 
branch points for a real system. In a straightforward ap- 
proach, one will have to locate those A inside the unit 
disk where Hx displays degeneracies. Although such a 
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program can be, at least in principle, carried out nu- 
merically, there are few chances of success without clues 
of where these points are located. This is because, in 
more than one dimension, these degeneracies occur, in 
general, at complex energies. One possible solution is to 
follow the lines presented in this paper: locate the type II 
degeneracies for a separable potential Vg, chosen as close 
to the real potential v as possible, and follow the trajec- 
tory of the branch points as the potential is adiabatically 
changed — Vs -\- j{v — w^), from 7 = to 1. We plan 
to complete such a program in the near future. 

The analytic structure of the band energies and Bloch 
functions of 3D crystals, viewed as functions of several 
variables kx , ky and kz is a much more complex problem, 
with qualitatively new aspects. It will be interesting to 
see if this problem can be tackled by the same analytic 
deformation technique. 

The main part of this work was completed while 
the author was visiting Department of Physics at UC 
Santa Barbara. This work was part of the "Near- 
sightedness" project, initiated and supervised by Prof. 
Walter Kohn and was supported by Grants No. NSF- 
DMR03-13980, NSF-DMR04-27188 and DOE-DE-FG02- 
04ER46130. The complex bands were investigated while 
the author was a fellow of the Princeton Center for Com- 
plex Materials. 



APPENDIX A 

We prove here that if {H\}\^c is an analytic family in 
the sense of Kator^ then i^m(A) defined in Eq. are an- 
alytic functions. If Rz{X) = (z — H\)^^, by definitionfi^ 
the limit 



lim 



RziX)-Rz{y) 
A- A' 



(Al) 



exist in the topology induced by the operator norm, for 
any A e C and z £ p{H\). We denote this limit by 
d\Rz{X). Consider now an arbitrary Aq G C, and a 
contour F in p{H\g) surrounding TV eigenvalues of -//aq- 
Note that the Bloch Hamiltonians considered in this pa- 
per have compact resolvent so their spectrum is always 
discrete. For A in a small neighborhood of Aq, F remains 
in p{H\) and we can define 



F^i\) = J^z"'Rz{\)^, 



2m 



(A2) 



and -F'm(A) = TrFm{X). Fm{X) is an analytic family of 
rank N operators for A in a small neighborhood of Aq. 
Indeed, if 



'dxRziX) 



dz 
2^' 



then 



F^{X) - Frn{X') 



A- A' 



(A3) 



(A4) 



as A' — > A, since it can be bounded by 



RziX)-Rz{X') 



A- A' 



dxRziX) 



\dz\ 
2n 



This means the limit 



lim 

A'^A 



F,n{X) — Fm{X') 

A- A' 



(A5) 



(A6) 



exists and is equal to (A) . Since (A) is the difference 
of rank N operators, it is at most rank 2A'^. In particular, 
|TrF;;,(A)| < 00. Then 



i^m(A)-F^(A') 



A- A' 







(A7) 



as A' ^ A, since 



Tr 



< m 



F^iX)-Fm{X') 
A - A' 

FraiX) - FraiX') 



A- A' 



Kix) 



(AS) 



and Eq. HA7|) follows from Eq. (jA4|) . Thus, the hmit 

hm (A9) 
A'^A A - A' 

exists and is equal to TrF'^-^{X). 

APPENDIX B 

We discuss here the analytic perturbations for linear 
molecular chains. As we did in the main text, we con- 
strain X and y in finite intervals. The Bloch functions are 
determined by the following Hamiltonian: 



H 



A,7 



-Aa + 7w, x,y £[(),b'] and z £[0,b], (Bl) 



where Aa is the Laplace operator with periodic bound- 
aries in X and y and the usual Bloch conditions in z. We 
show that if 



I'^IIl^ 



w(r)^c?r 



-,1/2 



< 00, 



(B2) 



then H\ .y is an analytic family for all 7 G C. For this we 
need the following technical result. 

Proposition. Suppose w satisfies Eq. (|B2|) . Then, for a 
positive and sufficiently large, there exists a positive Ca 
such that lim e„ = and: 



\\wf\\L^<ea\\{Hx,Q + a)f\\L-, 
for any / in the domain of H\xi. 



(B3) 
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Now, pick an arbitrary 70, let z € p{H\^^g) and de- 
note dz = ||(i?A,7o ~ ^)~^\\ < where |1 |1 denotes the 
operator norm. Since 

{I - eaho\)\\wf\\L2 < CaUH^^^, + a)fU2, (B4) 

for any / in the domain of H\,o, taking a sufficiently large 
so that 1 — ea|7o| > 0, we obtain: 



ea[l + \z + a\dz] 
1 - falTol 



< 00. 



(B5) 



If M denotes the right hand side of the above equation, 
then {H\,y ~ z)~^ is bounded for I7 — 70I < M^^ and has 
the following norm convergent expansion: 



A, 70 



zY 



(B6) 



Thus, (-ffA,7 ~ z) ^ is analytic at the arbitrarily chosen 
70- 

We now give the proof of the proposition. For / in the 
domainof Fa.o, let g = (i?A,o+a)/- If Ga = (i?A,o+a)"\ 
with a assumed sufficiently large so that {Hx q + a)^^ 
exists, we have 

/(r)= I Gx{v,v'-a)9{v')dv' (B7) 

and Schwartz inequality gives = sup 

r 

1/2 



< sup 

r 

If we denote 



j |GA(r,r';a)pdr' 



\9\W- (B8) 



aa = sup 



j \Gx{r,v'-a)\Hv' 



1/2 



(B9) 



with the aid of Eq. HB8|I , we obtain: 



\wf\\L-<aa\\w\\LA\{Hxfi + a)f\\L-, 



(BIO) 



i.e. Eq. IIB3|) . if we identify ea = oia\\w\\]^2. We remark 
that Ua defined in Eq. I|B9(I is optimal, in the sense that 
there are f's when we do have equality in Eq. ljB10|l . It 
remains to show that lima _> 00 aa = 0. 

If G" = (— A-|-a)~-^, with A the Laplace operator over 
the entire R'^, then we have the following representation: 

GA(r,r';a) = ^ A-"-G°(r - r' + R; a), (Bll) 
Rer 



where the sum goes over all points of the lattice F defined 
by R = {rixb' , Uyb' , Uzb). Since G*^ is real and positive, we 
can readily see that |GA(r,r';a)| < G|a| (r, r'; a). More- 
over, 

J G|A|(r,r';a)2dr'-(i/|A|,o + a)-'(r,r), (B12) 
and we have the following representation: 

(F|A|.o + a)-'(r, r) = ^ \X\--'C%R; a), (B13) 
Rer 

where G" = (—A -I- a)~^. Note that the diagonal part of 
{H\x\,o + is independent of r. G" can be exphcitly 
calculated, leading to: 



Rer 



j-v^|R| 



1/2 



(B14) 



with equality for A real and positive. The right hand side 
is finite for a sufficiently large and goes to zero as a — > 00. 
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